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SUMMARY 

^ The  oblique  entry  of  a two-dimensional,  fully  ventilated  foil 
into  a horizontal  layer  of  water  of  arbitrary  thickness  is  consid- 
ered. When  the  thickness  of  the  layer  is  finite,  exit  of  the  foil 
from  the  layer  is  studied.  The  present  work  is  an  extension  of  a 
previous  one^HTT^in  which  only  vertical  entry  and  exit  of  the 
foil  were  considered.  The  consideration  of  oblique  entry  makes  • 
the  present  solution  useful  in  studying  not  only  the  advance 
ratio  effect  on  a partially  submerged  supercavitating  propeller 
but  also  many  other  problems,  such  as  the  high-speed  water 
entry  and  exit  of  a wing  or  a control  surface,  the  operation  of  a 
high-speed  hydrofoil  in  a high  sea  state.  The  results  obtained  in 
this  work  indicate  that  the  force  and  moment  coefficients  of  a 
foil  decrease  as  the  advance  ratio  is  reduced. 
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1,  INTRODUCTION 

In  a previous  paper  [1] , the  present  author  used  a linearized 
flow  model  to  study  the  flow  characteristics  around  a partially 
submerged  supercavitating  propeller  operating  in  the  vicinity  of 
an  ocean  surface.  The  model  was  fashioned  after  the  blade  element 
theory  for  the  conventional  marine  propeller.  When  the  effects 
of  the  neighboring  blades  and  of  gravity  were  ignored,  the  flow 
resembled  that  of  the  entry  into  and  exit  from  a horizontal  layer 
of  water  of  finite  thickness  by  a two-dimensional,  fully  venti- 
lated foil  in  a gravity-free  space.  The  circular-arc  length  trav- 
eled in  the  ocean  by  the  leading  edge  of  the  blade  element  was 
taken  as  the  thickness  of  the  layer.  The  upper  and  lower  surfaces 
of  the  layer,  which  represented  the  ocean  surface,  were  assumed  to 
be  exposed  to  atmospheric  pressure.  The  foil,  which  represented 
the  blade  element,  entered  vertically  into  the  upper  free  surface 
of  the  layer  and  exited  from  its  lower  free  surface.  When  moving 
in  the  water,  the  foil  was  allowed  to  have  small,  unsteady  defor- 
mations, so  that  the  theory  could  be  applied  to  hydroelastic 
studies. 

In  studying  the  problem,  the  flow  motion  considered  in  [1J  was 
divided  into  three  different  phases:  the  initial  entry,  the  com- 
plete entry,  and  the  exit.  The  flow  motion  is  defined  as  being  in 
the  initial  entry  phase  when  the  foil  is  only  partially  submerged, 
in  the  complete  entry  phase  when  the  foil  is  completely  submerged, 
and  in  the  exit  phase  when  the  leading  edge  of  the  foil  pierces 
through  the  lower  free  surface  of  the  layer.  For  the  exit  phase, 
the  pressure  distribution  on  the  foil  was  explicitly  determined; 
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for  the  other  two  phases,  however,  it  was  determined  analytically 
except  for  a function  of  the  time  factor.  The  determination  of  the 
function  depended  on  the  solution  of  a system  of  integral  equations. 
The  integral  equations,  which  governed  the  flow  speed  on  the  cavity 
walls,  were  produced  as  a consequence,  when  the  role  played  by  the 
time  factor  in  solving  the  flow  velocity  was  reduced  from  an  ex- 
plicit variable  to  a parameter.  For  the  special  case  when  the 
thickness  of  the  layer  became  infinite  and  the  flow  was  in  its  ini- 
tial entry  phase,  solutions  for  arbitrary  foil  deformations  were 
obtained  ir  -ies  forms.  For  other  cases,  numerical  methods  were 
devised  to  solve  the  integral  equations. 

The  results  obtained  in  [1J  show  that  the  rate  of  increase  in 

the  loading  on  the  foil  is  most  rapid  during  the  initial  entry  and 

that  the  loading  almost  reaches  its  maximum  at  the  end  of  this 

phase.  During  the  same  phase,  the  influence  of  the  thickness  of 

the  layer  on  the  foil  loading  is  small  and  may  be  neglected  for 

most  cases.  This  finding  may  become  very  useful  when  the  struc- 
* 

tural  design  of  a propeller  blade  is  considered. 

The  theory  developed  in  (1],  however,  did  not  consider  the 
advance  ratio  effect,  which  is  known  to  have  some  influence  on  the 
thrust  and  torque  coefficients  of  the  partially  submerged  propel- 
ler. To  consider  the  advance  ratio  effect,  the  foil  should  enter 
the  horizontal  layer  of  water  obliquely,  instead  of  vertically  as 
in  tlj.  The  main  purpose  of  the  present  paper  is  to  consider  such 
an  effect:  that  is,  to  study  the  oblique  entry  and  exit  of  a fully 
ventilated  foil  from  a layer  of  water  of  finite  thickness.  The 
limiting  case  when  the  thickness  of  the  layer  becomes  infinite  is 
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also  considered. 

Numerical  results  for  a flat  plate  foil  and  for  circular-arc 
foils  are  obtained.  They  indicate  that,  when  the  angle  of  attack 
measured  from  the  entering  path  of  a foil  remains  constant,  the 
force  and  moment  coefficients  of  the  foil  decrease  as  the  advance 
ratio  is  reduced. 

2.  DESCRIPTION  OF  THE  PROBLEM 

When  the  blade  of  a partially  submerged  supercavitating  pro- 
peller enters  an  ocean  surface,  the  flow  is  expected  to  separate 
from  the  leading  edge  of  the  blade,  forming  on  the  suction  side  a 
cavity  that  vents  to  the  atmosphere.  When  the  flow  model  described 
in  [1]  is  adopted  and  the  advance  ratio  effect  considered,  the 
corresponding  flow  becomes  that  of  the  oblique  entry  of  a two- 
dimensional,  fully  ventilated  foil  into  a horizontal  layer  of  water 
of  finite  thickness,  as  shown  in  Fig.  1.  The  upper  and  lower  sur- 
faces of  the  layer  are  assumed  to  be  exposed  to  the  atmospheric 
pressure  p . The  foil  enters  the  upper  free  surface  and  exits  from 

a 

the  lower  one  with  a constant  velocity  U,  which  is  assumed  to  be 
large  enough  to  cause  the  flow  around  the  foil  to  become  separated 
but  not  so  large  as  to  render  the  incompressible  and  inviscid 
approximations  invalid.  In  this  work,  the  flow  separation  is 
assumed  to  occur  immediately  after  the  entry  of  the  leading  edge, 
and  the  pressure  inside  the  cavity  to  be  always  equal  to  p . 

a 

In  Fig.  1,  the  leading  edge  of  the  foil  is  denoted  by  B,  and 
the  trailing  edge  by  T.  The  cavity  walls  are  represented  by  the 
curves  BC  and  TA.  C and  A are  the  intersectional  points  where  the 
cavity  walls  and  the  upper  free  surface  of  the  layer  meet. 

When  the  chord  length  of  the  foil,  Jl,'  is  taken  as  the  charac- 
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teristic  length,  the  non-dimensional  thickness  of  the  layer  is 
denoted  by  h.  The  entering  speed  U and  time  £/U  are  taken  as  the 
characteristic  speed  and  time  respectively.  The  non-dimensional 
pressure  p is  defined  by 

P “ (p-Pa)/C>iPU2),  (1) 

where  p is  the  dimensional  pressure  and  p the  density  of  water. 
Based  on  these  characteristic,  quantities,  physical  parameters  of 
the  present  problem  can  be  made  into  non-dimensional  forms,  which 
are  used  in  the  following  analysis. 

For  reference,  a Cartesian  coordinate  system  (x,y)  is  chosen; 
the  x-axis  lies  on  the  original  undisturbed  upper  free  surface  of 
the  layer,  and  the  y-axis  points  vertically  upwards.  At  time  t^O, 
the  leading  edge  of  the  foil  is  assumed  to  touch  the  upper  free 
surface  at  the  origin  of  the  x-y  axes.  The  angle  between  the  ve- 
locity U and  the  x-axis  is  denoted  by  yir,  as  shown  in  Fig.  1. 
Without  loss  of  generality,  it  may  be  assumed  that  0<y<*j.  When 
the  advance  ratio,  X,  of  a bludf  element  is  defined  as  t'ie  ratio 
of  the  propeller's  advance  speed  to  its  circumferential  speed  at 
the  location  where  the  blade  element  is  situated,  X may  be  ex- 
pressed as 

X * tanpr y)ir.  (2) 

Suppor  during  its  motion  in  the  water  the  foil  is  performing 
small  (say,  of  the  order  of  e),  unsteady  motions  about  its  entering 
path,  and  suppose  the  percentage  change  in  the  chord  length  of  the 
foil  is  less  than  the  order  of  c.  Then,  to  the  first  order  in  e, 
the  locations  of  points  on  the  foil  surface  at  any  t may  be  de- 
scribed by  the  equation 


In  (3),  s is  the  distance  measured  from  the  origin  of  the  x-y  axes 
along  the  entering  path,  n the  distance  measured  normal  to  as 
indicated  in  Fig.  1,  and  g(s,t)  is  assumed  to  be  given,  e is  so 
small  that  the  conventional  linearized  approximation  used  in  solving 
an  inviscid  free  surface  flow  problem  for  small  disturbances  may 
be  used. 

Since  the  flow  is  considered  incompressible  and  inviscid, 
a velocity  potential  <j>(x,y,t)  may  be  used  to  describe  the  motion 
of  the  water.  Consequently,  when  the  effect  of  gravity  is  ne- 
glected, the  linearized  equation  of  motion  becomes 

d>t  + p/2  = 0.  (4) 

The  constant  pressure  condition  p=0  on  the  free  surface  may 
now  be  stated  as 

<J>t  - 0.  (5) 

Consistent  with  the  adopted  linearization  process,  this  already 
linearized  boundary  condition  need  not  be  applied  to  the  actual 
locations  of  the  free  surfaces.  Instead,  it  may  be  applied  to  the 
planes  y-0  and  y*-h,  and  to  the  portions  of  the  entering  path  that 
correspond  to  the  cavity  walls  as  e + 0. 

On  the  wetted  side  of  the  foil,  the  linearized  boundary  con- 
dition is 

♦n  s egt(s,t),  (6) 

which  is  applied  to  the  portion  of  the  entering  path  corresponding 
to  the  wetted  side  of  the  foil  as  e •*•().  At  x“t»,  which  are  denoted 
by  points  D and  E in  Fig.  1,  it  is  required  that 


1 


<J>X  = = 0 for  t>0.  (7) 

The  boundaries  to  which  the  linearized  boundary  conditions  are 
applied  define  the  linearized  flow  field.  Since  the  iayer  of  water 
is  assumed  to  be  originally  motionless,  the  initial  condition  may 


be  stated  as 


in  the  entire  flow  field. 


<p  = 0 


As  explained  in  HI,  it  is  possible  to  express  the  boundary 


condition  (5)  as 


* 0 


on  the  free  surfaces  y=0  and  y=-h,  and  as 

<t>s  = <l(s) 


(10) 


on  the  cavity  walls  where  n=0.  In  (10),  q(s),  the  speed  function 
on  the  cavity  walls,  is  not  yet  known,  but  it  will  be  determined 
in  the  process  of  solving  the  problem.  The  advantage  of  replacing 
the  free  surfac  boundary  condition  (5)  by  (9)  and  (10)  is  that 
the  variable  t m acts  merely  as  a parameter  in  solving  the  flow 
velocity. 

3.  SOLUTIONS 

The  boundary  conditions  (6),  (7),  (9)  and  (10)  form  a typical 
boundary-value  problem,  whose  solution  determines  the  velocity 
field  that  includes  the  unknown  speed  function  q(s).  To  solve  the 
present  problem,  the  same  methed  of  solution  used  in  [1]  is  adopted 
here.  The  flow  motion  is  divided  into  three  different  phases:  the 
initial  entry,  the  complete  entry,  and  the  exit.  The  solution  to 
each  phase  is  sought  separately. 
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To  find  the  solutions  for  the  initial  and  the  complete  entry 
phases,  the  linearized  flow  field  is  mapped  onto  the  upper-half 
of  a parametric  c -plane  for  a given  t,  as  shown  in  Fig.  2,  where 
the  images  of  points  A,  B,  C,  D,  E and  T are  labelled  accordingly: 
point  A is  mapped  to  the  origin  of  the  £-plane,  B to  (b,0),  C to 
(c,0),  D to  the  point  of  infinity,  E to  (-1,0)  and  T to  (t,0). 

The  mapping  function  may  be  expressed  as 


ff  - c‘Y(e-b)U-c)Y'1(c*i)"1, 


where  z=x<-iy  and  c-5+i n.  In  (11),  the-  branch  cut  is  chosen  to  lie 
between  £=0  and  £=c  with  0 £ Arg(c)  < 2it  and  0£Arg(?-c)  < 2tt.  It 
can  be  shown  from  (11)  that 

1 

c = (l+b)^”^-l  (12) 

and  that  b as  a function  of  t is  governed  by  the  equation, 

rt/h  = /VY(l-y)(c/b-li)Y'1(l/b+u)‘1du.  (13) 

0 

The  integral  shown  in  (13)  can  be  expressed  in  terms  of  a hyper- 
geometric function  of  two  variables;  however,  it  is  more  convenient 
for  future  applications  to  express  the  relationship  between  b and 
t numerically.  The  results  for  several  values  of  X are  shown  graph- 
ically in  Fig.  3.  It  should  be  noted  that  the  correspondence  of  b*0 
with  t=0  and  of  b*»  with  t*hcscyir  can  each  be  shown  analytically. 

On  the  boundary  0£££c,  the  dependence  of  s on  £ and  t may 
be  expressed  in  an  integral  form  as 

A -V V-  1 . 


irs/h  « / v”Y(b-v)  (c-v)Y’*(l+v)  ’*dv. 
0 


(14) 


Obviously,  s has  only  e~e  extremum  in  the  region  0<££c.  This 
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occurs  at  £=b,  at  which  s=t.  As  £ increases  from  0 to  b,  s in- 
creases monotonically  from  0 to  t;  as  5 further  increases  from  b 
to  c,  s decreases  monotonically  from  t to  0. 

Initial  Entry 

This  phase  of  motion  occurs  during  the  time  interval  0£t_<l, 
in  which  the  trailing  edge  of  the  foil  has  not  yet  become  submerged. 
Consequently,  only  one  cavity  wall  is  generated.  This  is  located 
on  the  suction  side  of  the  foil  and  is  described  in  Fig.  1 as  BC. 

Suppose  the  unknown  speed  function  q(s)  on  the  cavity  wall 
BC  is  designated  by  q^(s).  By  solving  the  boundary-value  problem 
(cf.  fl]),  the  complex  velocity  w(z,t),  which  satisfies  the  finite 
pressure  condition  at  point  A and  C and  the  %-root  type  of  singu- 
larity in  pressure  at  point  B,  becomes 

w * iT  -fcc)1  \ (At(s,t)*/Cqi(s)|^-%^  AL'  (IS) 

"l  C Y(c-b)*  0 1 b 1 (c-{)  Y e-C 

where  q^(s)  is  the  solution  of  the  integral  equation 


t r % 

/ Igt  (s , t)+q,  (s)  J C — - — ) ds  « 0.  (16) 

0 1 |£-b| 

In  the  above  expressions,  C is  regarded  as  a function  of  z and  t, 

£ a function  of  s and  t,  and  the  function  (c-b)  is  restricted  to 
the  branch  specified  by  0£Arg(^-b)  < 2ir.  It  should  be  noted  that 
e,  the  parameter  used  to  indicate  the  smallness  of  the  disturbances, 
has  been  omitted  from  (15)  and  (16).  The  omission  of  e will  con- 
tinue from  here  onward. 

When  gt(s,t)  is  given,  the  integral  equation  (16)  determines 
q^s)  for  0<s<t£l.  However,  for  the  limiting  case  as  t + 0, 
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(17) 


q^(0)  can  be  explicitly  determined,  and  the  result  is 

qx  C°)  = -Q0(y)gt(o,o), 


where 


and 


1 /T^y 

Qa(y)  = / K(y)dp//  K(p)dp, 
U 0 / 1 


K(u)  = M%-Y|l.-M|li(Ii7-M)r'1. 


(18) 


(19) 


Equ.  (17)  indicates  that  q ^ (0 ) is  independent  of  h,  the  thickness 
of  the  layer.  The  function  Qq(y)  is  plotted  in  Fig.  4,  where  the 
ordinate  is  Qq  and  the  abscissa  y>  whose  corresponding  value  of  A 
is  also  shown  in  the  figure.  With  the  aid  of  either  (9)  or  (15), 
it  can  be  shown  that  the  velocity  at  point  C is  always  pointing 
vertically  upward,  irrespective  of  the  value  of  y so  long  as  0 < yf. 
h.  The  magnitude  of  this  velocity  is  equal  to  Qq (Y)gt (0,0)cscyn . 

Throughout  the  initial  entry  phase,  A (i.e.,  point  c-0)  is 
the  point  of  contact  between  the  foil  and  the  upper  free  surface 
of  the  layer.  The  asymptotic  behavior  of  w in  the  neighborhood  of 
A can  be  studied  from  (15).  For  the  case  when  0 < y < Hf 

w 'v  a^(t)cY_5*  as  C^O,  (20) 

where 

a,(t)  - ^e'iY*{/bgt(s,t)*/Cq1(sm'’’'1,K-b|'‘(c-0Y'1d5.C21) 

1 ib’  0 1 b 1 

When  ysh,  however, 

w * |gt(0,t)  In  c as  C + 0.  (22) 

The  pressure  distribution  on  the  wetted  side  of  the  foil, 
yhich  lies  on  the  segment  0£C<.b,  can  be  derived  from  (15).  The 
derivation  involves  the  integration  of  w(z,t)  with  respect  to  z 
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(which  produces  the  complex  velocity  potential  f(z,t))  and  the 
differentiation  of  f(z,t)  with  respect  to  t.  The  final  result 
may  be  written  as 


ln 


[✓cTT-TWg(b-g)}' 

bU-C| 


ds , 


(23) 


where 


a2(t)  = {/*t(s,t)-/Cq1(s)}e%"1r|5-b| 


'!i(c-C)1,'1d5. 


(24) 


The  above  results  form  the  basis  on  which  numerical  calculations 
for  the  solution  of  the  initial  entry  can  be  made. 

Complete  Entry 

This  phase  occurs  during  the  time  interval  lctchcscyn,  in 
which  both  the  leading  and  trailing  edges  of  the  foil  remain  sub- 
merged in  the  water.  In  this  case,  a second  cavity  wall,  TA,  is 
generated  in  addition  to  the  already  established  BC.  The  flow 
configuration  has  been  shown  in  Fig.  1.  Suppose  the  unknown 
speed  function  on  BC  where  s > 1 is  again  designated  as  q^fs)  and 
that  on  TA  as  q2(s).  Then,  the  solution  for  w(z,t)  may  be  written 
as  (cf.  Ill) 

i r-c  1_Y  r-T  **  T 

b C r 1-Yir  wl1*  Ac 

♦/  gt(s,t)*/  qi(s)>(^)  (25) 

where  t is  the  ^-coordinate  of  the  trailing  edge  T,  whose  s -coor- 
dinate is  (t-1),  and  the  function  (c-t)1*  is  restricted  to  the 
branch  specified  by  0 < Arg(c-t)  < 2w.  The  value  of  t can  be 
determined  from  (14)  by  iterations.  The  iteration  process  is  a 


straightforward  one  and  can  be  made  to  converge  for  any  s on  the 
foil  or  on  the  cavity  walls.  This  is  because  s is  a monotonically 
increasing,  or  decreasing,  function  of  £ in  those  ranges. 

The  unknown  speed  functions  q^s)  for  1 < s < hcscyir  and  q2(s) 
for  0 < s < hcscyn- 1 are  governed  by  the  following  two  integral  equ- 
ations , 


t-1  t t , 

{/  q7(s)+  / gt(s,t)+/  q-.  (s) }- ~j~ 

0 2 t-1  1 0 1 I (5-t) (C-b) I’ 


1 

0 

ds  * 

£ 

0 

(26) 


Using  the  procedures  described  in  the  previous  section,  the 
pressure  p on  the  foil  BT,  which  lies  on  the  segment  x<  c_<b,  can 
be  obtained  from  (25),  and  the  result  may  be  expressed  as 


P ■ 


2 r-b  1_Y  r-r  ** 


| /Vrs.t) i, sm-u snumm rEsnlis_ 

11  * - tt:  (b-x)|c-5| 


t-1 


where 

a3(t) 


{/ 

0 


1-Y 


d£ 


q;.(s )+/  gt(s,t)-/  q1(s)}(rrr)  1 7“ 75 

" x t b 1 (5-b)(C-x)r 


(27) 


. (28) 


When  the  shape  function  g(s,t)  is  given,  (25),  (26)  and  (27)  can 
be  used  to  calculate  the  solution  of  the  complete  entry. 

Special  Case  as 

As  h-*“,  the  problem  becomes  that  of  the  oblique  entry  of  a 
fully  ventilated  foil  into  a deep  ocean.  The  solutions  for  the 
present  case  may  be  extracted  from  those  obtained  previously  by 
re- scaling  the  parametric  5 -plane  in  such  a way  that  point  E tends 
to  point  D at  the  point  of  infinity.  This  can  be  achieved  by 
replacing  b,  c and  x in  the  previous  solutions  by  c/h,  b/h,  c/h 
and  x/h,  respectively,  and  then  letting  h *► 09 . When  such  a process 
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is  carried  out,  (11)  becomes 


H - -i?'YU-b)U-c)Y'1, 


and  (12)  and  (13)  become,  respectively, 

c = b/(l-y) 


where 


Similarly,  (14)  becomes 


b = B(y)irt, 


1-y  Y 

b(y)  - (y)  • 


where 


r 

s * — / v'Y(b-v)(c~v)Y'1dv  = tS(y), 
n 0 


s Cu)  3 Bfy)/^^-  m)y 


y 3 £/b.  (35) 

For  both  the  initial  and  complete  entries,  the  solutions 
w(z,t),  the  integral  equations  which  govern  qj(s)  and  q2(s)>  and 
the  pressure  distributions  on  the  foil  retain  formally  the  same 
expressions  as  those  shown  in  the  previous  sections  where  the  case 
of  finite  h was  studied,  with  the  exception  that  the  relations  (11) 
to  (14)  should  be  replaced  by  (29),  (30),  (31)  and  (33). 

A special  case  for  which  a closed  form  solution  can  be  ob- 

4 

tained  during  the  initial  entry  is  when  gt(s,t)  is  given  by 

gt(s,t)  » t0sa,  (36) 

where  $ and  a are  arbitrary  constants.  To  seek  the  closed  form 
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solution,  the  integral  equation  (16)  is  transformed,  by  using  (30), 
(33)  and  (35),  into 

1 

1 T-y 

{/  gf (s,t)'+/  q1(s)}K(y)du  = 0,  (37) 

0 z 1 1 

where  K(y)  is  given  by  (19).  The  substitution  of  (33)  into  (36) 


gives 


gt(s,t)  = te+aSCT(u) . 


When  (38)  is  substituted  into  (37),  it  becomes  obvious  that  the 
solution  of  the  integral  equation  (37)  is 

„ _ a „&+0  rxQ' 


qi  (s) 


Vs 


Where  A„  is  a constant, 

pO 


1 

Aq_  * /1Sa(u)K(u)dy//r^S3+cf(u)K(M)du. 


Note  that  Agg^C^y).  When  the  motion  of  a foil  can  be  expressed 
as  a linear  combination  of  terms  of  the  type  shown  in  (36),  the 
initial  entry  solution  is  readiiy  obtainable  by  the  aid  of  (39). 
Examples  have  been  given  in  [1],  in  which  the  heaving  and  pitching 
motions  of  a vertically  entering  flat  plate  were  studied. 

After  complete  entry  is  achieved,  the  foil  is  assumed  to 
move  indefinitely  into  the  deep  ocean  without  changing  its  phase 
of  motion.  Eventually,  as  t + »,  the  flow  tends  to  the  classical 
cavity  flow  having  a zero  cavitation  number,  and  no  exit  phase  will 
occur  for  this  special  case.  For  the  case  of  finite  h,  the  exit 
phase  is  considered  in  the  following  section. 


EXIT 


As  the  leading  edge  of  the  foil  pierces  through  the  lower  free 
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surface,  the  layer  of  water  is  assumed  to  be  separated  into  two 
disconnected  flow  regions  bounded  by  ATB'EA  and  DB"CD,  respectively, 
as  shown  in  Fig.  5.  In  the  figure,  B*  indicates  the  point  of  con- 
tact between  the  foil  and  the  lower  free  surface,  and  B”  the  point 
at  which  the  cavity  wall  and  the  lower  free  surface  intersect. 

All  other  points  in  the  figure  have  previously  been  defined. 

To  study  the  hydrodynamic  forces  acting  on  the  foil  during  the 
exit  phase,  which  occurs  when  hcscyTr <_  t < hcscyrr+l , only  the  flow 
region  bounded  by  ATB'EA  needs  to  be  considered.  It  should  be  men- 
tioned here  that  the  entire  flow  motion  from  t=0  to  t=hcscyTi+l  is 
continuous,  since  the  flow  condition  at  the  end  of  each  phase  of 
motion  is  used  as  the  initial  condition  for  the  next  phase  of  motion. 

When  the  linearized  approximation  is  taken,  the  linearized 
flow  field  can  be  mapped  onto  the  upper-half  of  a parametric  C-  , 
plane,  as  shown  in  Fig.  6,  with  point  A mapped  to  the  origin  of  the 
C-plane,  point  E to  1 and  point  B*  to  the  point  of  infinity. 

The  mapping  function  can  be  extracted  from  (11)  by  letting  b-*-®, 
and  the  result  is 

|f-  -le^V^ir1.  (41) 

The  boundary  conditions  for  the  present  phase  of  motion  may 
be  stated  as  follows: 


- 0 

on 

ti 

o 

+ 

>• 

K<  0, 

(42) 

q2^  35 

on 

23 

II 

O 

+ 

o<  e < t , 

(43) 

■ M8’*5  lf| 

on 

na0+, 

r<  e 

(44) 

A _ = <b  » 0 

yn 


at  5=  - 1 , 


where  s as  a function  of  £ may  be  expressed  as 

s - ~ /V^S+D^di-  (46) 

* 0 

and  t,  again,  is  the  ^-coordinate  of  the  trailing  edge  T.  As  men- 
tioned before,  the  flow  condition  at  the  end  of  each  phase  of 
motion  should  be  used  as  the  initial  condition  for  the  next  phase. 
This  result  is  satisfied  if  q2(s)  in  (43),  when  s<hcscy7r-l,  is  the 
solution  found  at  the  end  of  the  complete  entry. 

The  boundary-value  problem  specified  by  (42)  to  (45)  can  be 
solved,  and  the  solution  which  satisfies  the  finite  pressure  con- 
dition at  points  A,  T and  B'  is 

w - ie’1'|,,ri;Y(C-T)ll{/Tq2(s)*/gt(s,t)}5'Yn-Tr)i  0^  , (47)- 

subject  to  the  condition  that 
t-1  hcscyv 

{/  q7(s)+  / g (s,t)}|5-T|“^ds  « 0,  • (48) 

0 1 t-1  z 

where  0<  Arg(5-r)  < 2v  has  been  chosen. 

In  (48),  q2(s)  is  considered  as  known  when  0 £ s < hcscyTT-1 , 
and  therefore  (48)  is  an  integral  equation  for  q2(s)  only  when 
hcscyTt - 1 £ s < hcscyir.  When  this  part  of  q2(s)  becomes  known,  w(z,t) 
for  the  exit  phase  can  be  determined  from  (47).  However,  the  as- 
ymptotic behavior  of  w at  point  B',  where  Id"®0,  can  be  visualized 
directly  from  (47)  without  solving  (48);  the  result  is  that  w at  B’ 
is  bounded  if  y < H and  becomes  logarithmically  singular  if  y=*s. 

An  interesting  aspect  of  flow  behavior  is  revealed  by  studying 
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isssfe 


the  velocity  potential  f(z,t).  When  w(z,t)  given  in  (47)  is  inte- 
grated with  respect  to  z,  f(.£ , t)  becomes 

2 


i t-1  hcscyir  vfTi  * 

f“fT'?{ig  q2^s^+t,/’1  g t ( s , t ) } 1 n-[-  .5 ' i " J—  d s , 


(49) 


where  fT  is  the  value  of  f(z,t)  evaluated  at  the  trailing  edge  T, 
5 is  regarded  as  a function  of  s and  t,  and  c as  a functirr  of  z 
and  t.  The  partial  differentiation  of  f with  respect  to  t by 
holding  z fixed  yields 


8f 


-I  / 

8t  *t-l 


hcscyit 


2 cl  i|i 

gtt(s,t)lni^|i^Il_ds  gt(t-l,t)l  , (50) 


where  i|>T  is  the  stream  function  evaluated  at  T.  From  (50),  it  is 
seen  that  <j>t=0  identically  if  gtt=0.  This  happens,  for  example, 
when  the  foil  is  a flat  plate.  In  such  a case,  the  potential  <|> 
becomes  stationary,  and  therefore  the  pressure  becomes  zero  through- 
out the  entire  flow  field.  In  other  words,  for  a foil  whose  gtt**0 
during  exit,  the  flow  around  the  foil  becomes  conformed  to  the 
motion  of  the  foil  so  that  no  additional  disturbance  is  generated. 

To  exit,  the  foi1  simply  slides  out  of  the  already  established  flov 
without  encountering  any  resistance. 

For  a general  case  when  gtt(s,t)^0,  the  pressure  acting  on 

the  wetted  part  of  the  foil,  which  is  represented  by  TB*  where 
5 > t , is 


| /hcscY\t(s,t)lnMHI-EIlIds. 

* t-i  " U-?| 


(51) 


The  determination  of  p does  not  depend  on  the  solution  of  the  inte- 
gral equation  (48).  To  determine  w(z,t),  however,  (48)  has  to  be 
solved. 
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3.  EXAMPLES 

Consider  a circular-arc  foil  with  a small  arc  angle  20  and  a 
small  angle  of  attack  a from  the  entering  path.  In  this  case, 

g(s,t)  = (a-0)(t-s)+0(t-s)2.  (52) 

This  foil  degenerates  into  a flat  plate  as  0 -*•().  When  (52)  is 
substituted  into  (16)  and  (26),  the  unknown  speed  functions  q^(s) 
and  q2(s)  can  be  solved  from  the  resulting  integral  equations. 

The  method  of  solution  has  been  explained  in  detail  in  [1] ; for 
convenience,  however,  it  is  briefly  described  in  the  following. 

Since  q^(s)  is  expected  to  be  bounded  everywhere  on  the  cavity 
wall  BC,  it  may  be  expanded  into  a regular  series  as 

CO 

q.(s)  = l b (s-s  )n  (53) 

* n*0  mn  m 

for  s in  the  neighborhood  of  s=sm  (m=0,  1,  2,  ...),  which  is  a set 
of  fixed  points  on  the  cavity  wall  BC.  The  bmn  are  unknown  con- 
stants to  be  determined. 

A regular  expansion  like  (53),  however,  cannot  be  made  for 
q2(s)  because  it  is  singular  at  s=0.  But,  when  the  singular  part 
is  removed,  the  remainder  may  be  expanded  into  a regular  series. 
Suppose  the  singular  part  of  q,(s)  is  denoted  by  q|(s).  Then,  from 
(20)  and  (14),  q£(s)  may  be  expressed  as 

q|(s)  - -C(l)s”(,s"Y)/(1‘Y)  when  y < h,  (54) 


where 


crt1  . 1.  h .(,1‘T)/(1‘r)cV>>Y/(l-Y) 

C(t)  c 6 

b c 

•{/  gt(s,tw  q1(s)H'!,'YIC-b|’,(c-5)r'1«. 
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qf(s)  3 ^gt(°»l)ln when  y=h.  (56) 

Equ.  (56)  is  purposely  constructed  in  such  a way  that  q£(s)  tends 
to  zero  as  s-+°>.  It  can,  therefore,  be  applied  even  to  the 
special  case  when  h ■+  Because  the  trailing  edge  of  the  foil  is 
always  one  unit  above  the  leading  edge,  [q2 (s) -q| (s) ] may  conven- 
iently be  expanded  into 

q2Cs)-q*(s)  = J°  cm„Cs-yD\  (57) 

where  c are  unknown  constants, 
mn 

To  study  the  initial  entry  phase,  (53)  is  truncated  to  a poly- 
nomial of  degree  N and  substituted  in  (16).  The  integral  equa- 
tion is  now  ready  to  be  solved  in  steps.  Suppose,  when  m=»0,  sQ=tQ 
=0  are  chosen.  Then,  (16)  contains  bQn  (n=0,  1,  2,  N)  as  un- 

known constants.  These  unknowns  can  be  determined  by  evaluating 
(16)  at  N+l  points  in  t,  starting  from  the  first  point  at  t=tQ+6, 
where  6 is  a small  number,  and  increasing  consecutively  to  the 
(N+l)  point  denoted  by  t*tp  which  is  assumed  to  be  less  than  1. 

At  these  points  in  t the  integral  equation  is  satisfied.  The 
resulting  N+l  algebraic  equations  determine  the  bgn,  and  therefore 
q^(s)  becomes  known  for  0<_  s<_  tp  Since  q^(s)  is  now  known  when 
(Hs^tp  the  portion  of  the  integral  equation  associated  with 
it  may  be  transfered  to  the  right-hand  side  of  tha  equation  so  that 
the  foregoing  process  may  again  be  used  to  obtain  the  solution 
qj (s)  when  t ^ < s < 1. 

Calculations  were  carried  out  for  a flat  plate  foil  at  an 
angle  of  attack  a and  for  a circular-arc  foil  with  0»o.  The  ad- 
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vance  ratios  chosen  were  A=0,  0.5  and  1.  In  the  calculations,  the 
relationship  between  b and  t was  determined  ro;i  (13),  which  has 
been  shown  graphically  in  Fig.  3.  The  dependence  of  s on  £ and  t 
was  determined  directly  from  the  integral  (14).  But,  for  the  in- 
verse case  when  s and  t ..ere  prescribed,  the  corresponding  £ was 
determined  from  (li)  by  iterations.  As  explained  before,  the  iter- 
ation process  is  a straightforward  one  and  can  be  made  to  converge 
for  any  s on  the  foil  or  on  the  cavity  walls.  For  the  initial 
entry  phase,  the  value  of  h need  not  be  specified  in  the  calcula- 
tions, and  it  is  considered  a free  parameter  of  the  problem. 

To  study  the  complete  entry,  which  occurs  during  the  time 
interval  l<t<hcscyir,  the  integral  equations  (26)  --  instead  of 
(16)  --  should  be  solved.  During  this  phase  of  motion,  q2(s)  is 
an  unknown  function  in  addition  to  q^s)  when  s > 1.  When  the 
series  expansion  (57)  for  [q2(s)-qJ(s)J  is  truncated  to  a poly- 
nomial of  degree  N,  the  foregoing  method  of  solution  may  again  be 
used  to  determine  q^s)  when  1 < s < hcscyir  and  to  determine  q2(s) 
when  0 < s < hcscyir- 1.  In  the  present  case,  h cannot  be  considered 
a free  parameter  as  it  was  during  the  initial  entry  because  the 
value  of  h is  needed  in  the  determination  of  t,  the  ^-coordinate 
of  the  trailing  edge  of  the  foil. 

The  pressure  distributions  on  the  foil  during  the  initial  and 
complete  entries  can  be  determined  from  (23)  and  (27),  respectively. 
Their  properties  are  similar  to  those  shown  in  [11 , where  the  spe- 
cial case  of  y*h  was  studied.  For  illustration,  some  of  the  pres- 
sure distributions  are  shown  in  Fig.  10.  The  explanation  of  the 
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figure  will  be  given  at  a later  stage  in  this  work. 

As  explained  previously,  closed  form  solutions  can  be  obtained 


for  the  initial  entry  phase  as  h + 00 . When  (52)  is  compared  with 


(36),  the  solution  q^s)  can  easily  be  obtained  by  the  aid  of  (39) 


The  result  is 


qx(s)  =*  - (a-0)Qo(Y)-20Q1(y)s, 


where  Qj (y)=A10-A01  and  is  defined  in  (40).  The  pressure 


distribution  on  the  foil  is 


p = [(a-0)Po(y)+0tP1(y)  ] [p/(l-p)K  +0tG(p)  for  0<p<l.  (59) 


Where  p is  defined  in  (35), 


po«  ■ 


1 _i_ 

pi(Y)  3 tl-S(y)]  +Q1(y)/177S(u)}-^M-  dp  (61) 

1 ni  YJ  0 1 1 1-p 


G(p)  = f ( 1 - v ) / v ] (v) In——  (1~ /v 1 1 BUll  dv . 

n 0 ( p~v | 


The  function  G(p)  is  always  positive  when  0<p<l  and  vanishes 


when  ps0  or  p*l.  Since  0£S(p)£l  when  0£p£l,  it  is  not  diffi- 
cult to  show  that  Pq(y)  and  P^(y)  are  positive  when  0<y<_*j.  This 
implies  that,  if  0<6_<al  the  pressure  near  the  leading  edge  B of 
the  foil  is  always  positive  definite  during  the  entire  initial  entry 


phase,  a condition  that  will  prevent  ventilation  from  developing 


on  the  pressure  side  of  the  foil.  For  reference,  the  function 


Ql(y),  Pq(y)  and  P^(y)  are  shown  graphically  in  Fig.  7,  where  the 
scale  for  is  shown  on  the  left-hand  side  of  the  figure  and  that 
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for  PQ  and  on  the  right-hand  side. 

When  0=0,  the  foregoing  results  correspond  to  those  of  a flat 
plate  foil.  In  this  case,  p shown  in  (59)  is  a function  of  y only, 
i.e.,  a function  of  the  combined  variable  s/t  only.  This  means 
that,  as  expected,  a similarity  solution  exists  for  a flat  plate 
during  the  initial  entry  when  h=». 

For  the  complete  entry,  however,  no  closed  form  solution  can 
easily  be  obtained  even  for  the  case  of  h=«.  The  numerical  method 


described  previously  was  used  to  obtain  the  solution. 

For  a finite  h,  exit  occurs  during  the  time  interval  hcscyv £ 
t<hcscyir+l.  For  this  phase  of  motion,  the  pressure  distribution 
on  the  foil  can  be  determined  directly  from  (51).  In  the  calcula- 
tions, the  functional  relationship  between  s and  £ was  obtained 
from  (46):  s was  calculated  directly  from  (46)  when  5 was  given, 
but  £ was  calculated  by  iterations  when  s was  given.  Some  of  the 
pressure  distributions  during  the  exit  phase  are  shown  in  Fig.  10. 

The  force  and  moment  coefficients  of  the  foil  are  obtained 
by  integrating  p,  which  is  defined  in  (1),  over  the  wetted  part  of 
the  foil.  Suppose  the  force  coefficient  in  the  n-direction,  when 
normalized  by  the  chord  length  of  the  foil,  is  denoted  by  Cn>  Then, 

Cn  - -/pds.  (63) 

The  moment  coefficient  about  the  leading  edge  of  the  foil,  CM  , 


may  be  expressed  as 


' * /(t-s)pds, 

mB 


where  CM  is  positive  when  the  moment  is  counterclockwise.  For 
the  particular  type  of  foils  under  consideration,  the  force  coef- 
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(65) 


ficicnt  in  the  s-direction,  C,.,  is  related  to  Cn  and  CM  by 


JM 


B 


Ce  " (a - 0 ) Crt  ■ 20Cj^  . 


n 


B 


Numerical  results  of  C and  Cw  for  a flat  plate  foil  at  an 

n mb 

angle  of  attack  a are  shown  in  Fig.  8.  Fig.  8(a)  shows  the  results 
of  *Cn/a  versus  t for  various  values  of  h.  For  each  value  of  h, 
a group  of  three  curves  is  presented  in  Fig.  8(a);  they  correspond 
to  different  values  of  X,  as  indicated  in  the  figure.  Generally 
speaking,  when  other  factors  are  equal,  the  loading  is  larger  when 
h is  greater,  and,  within  a group  of  constant  h,  the  loading  is 
larger  when  X is  larger.  Fig.  8(b)  shows  the  corresponding  curves 
for  . For  comparison,  the  force  and  moment  coefficients  of  the 
corresponding  supercavitating  plate,  at  zero  cavitation  number, 
are  also  shown  in  Fig.  8(a)  and  (b)  at  the  location  indicated  by 
t=®.  They  are  the  upper  bounds  of  the  respective  coefficients  in 
the  present  study. 

Since  the  force  (or  moment)  coefficients  for  different  values 
of  h and  X are  nearly  equal  during  the  initial  entry  (0<t£l), 
they  cannot  be  distinguished  clearly  in  Fig.  8(a)  (or  8(b)).  For 
a clearer  view,  they  are  plotted  in  different  scales  in  Fig.  8(c). 
In  reading  Fig.  8(c),  t should  be  considered  less  than,  or  equal 
to,  unity  for  any  given  set  of  values  of  h and  X. 

Similar  results  for  a circular* arc  foil  with  8«a  are  shown  in 
Fig.  9. 

Some  representative  pressure  distributions  on  the  foils  are 
shown  in  Fig.  10.  Fig.  10(a)  shows  those  on  the  flat  plate  foil. 
Fig.  10(b)  on  the  circular-arc  foil,  and  Fig.  10(c)  on  the  circu- 
lar-arc foil  during  the  exit  phase. 
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The  four  solid  curves  shown  in  Fig.  10(a)  indicate  the  pres- 
sure distributions  for  the  case  when  h=°°  at  various  times  and  dif- 

'K 

ferent  values  of  X.  When  the  time  factor  alone  is  considered, 
these  curves  may  be  divided  into  two  groups : one  corresponds  to 
t=«>,  and  the  other  to  0<t£l,  as  indicated  in  the  figure.  The 
latter  group  represents  the  pressure  distributions  on  the  flat 
plate  during  the  entire  initial  entry  phase  (0<t£l)  for  different 
values'  of  X as  indicated.  The  loading  is  larger  when  X is  larger. 
The  former  group,  which  consists  of  a single  solid  curve,  repre- 
sents the  pressure  distribution  on  the  plate  at  any  X when  t=«. 

This  pressure  distribution  is  actually  the  same  as  that  on  the 
corresponding  supercavitating  plate  at  zero  cavitation  number. 

These  two  groups  of  curves  form  the  upper  and  lower  bounds 
for  the  pressure  distributions  on  a flat  plate  foil  entering  an 
infinitely  deep  ocean  at  one  of  the  values  of  X considered.  In 
fact,  the  former  group  is  the  upper  bound  for  all  possible  pres- 
sure distributions  obtainable  in  tK;  present  study.  In  conjunc- 
tion with  Fig.  8(a),  these  two  grow  s are  sufficient  for  visual- 
izing the  pressure  distributions  on  & flat  plate  entering  a layer 
of  water  of  finite  thickness. 

Suppose,  for  example,  the  case  of  ha2ir  is  considered.  During 
the  initial  entry,  the  loading  on  the  plate,  as  shown  in  Fig.  8(a), 
is  about  the  same  as  that  for  the  case  of  h®«.  This  means  that, 
for  those  values  of  X considered,  the  pressure  distributions  on 
the  plate  during  the  initial  entry  are  about  the  same  as  those 
shown  in  the  group  indicated  by  0<t<_l  in  Fig.  10(a).  During  the 
complete  entry,  the  loading  on  the  plate,  as  shown  in  Fig.  8(a), 
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increases  continuously  from  that  at  t=l  until  t*3.1,  at  which  time 
the  maximum  loading  on  the  plate  is  about  reached.  This  indicates 
that  during  the  time  interval  l£t£3.1  the  pressure  distributions 
on  the  plate  shift  upward  from  one  of  the  group  indicated  by  0<t£ 
1 shown  in  Fig.  10(a)  to  the  corresponding  one  in  the  group  of 
three  broken  curves,  shown  in  the  same  figure,  according  to  the 
value  of  A considered.  After  t=3.1,  the  loading  on  the  plate 
starts  to  decrease.  Therefore,  the  pressure  distributions  on  the 
plate  start  to  shift  downward  from  the  group  of  three  broken  curves 
shown  in  Fig.  10(a),  and  eventually  they  become  identically  zero 
during  exit. 

For  the  case  when  h=0.5n,  the  situation  is  slightly  different. 
From  Fig.  3(a),  it  is  seen  that  the  loading  on  the  plate  has  about 
reached  its  maximum  at  the  end  of  the  initial  entry  phase  (t*l) . 

But  this  maximum  loading  is  less  than  the  loading  for  the  case  of 
h=»  at  t*l.  This  means  that  for  the  case  of  h=0.5ir  the  pressure 
distributions  on  the  plate  for  the  values  of  A considered  never 
exceed  those  shown  in  the  group  indicated  by  0<t£l  in  Fig.  10(a). 
When  t is  small,  the  pressure  distributions  are  about  the  same  as 
those  shown  in  the  group.  As  t becomes  larger,  but  still  less  than 
unity,  the  pressure  distributions  over  a major  part  of  the  plate 
start  to  shift  downward  from  those  shown  in  the  group  except  near 
the  trailing  edge  of  the  plate,  where  the  pressures  may  actually 
shift  upward  for  a while,  and  eventually  all  the  pressure  distri- 
butions become  identically  zero  during  exit. 

For  the  circular-arc  foil,  similar  situations  occur,  although 
the  pressure  distributions  during  exit  are  not  identically  zero. 
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Some  typical  pressure  distributions  are  shovm  in  Fig.  10(b)  and 
(c) . Fig.  10(b)  is  equivalent  to  Fig.  10(a).  Fig.  10(c)  shows 
the  pressure  distributions  on  the  foil  during  exit  for  the  case 
of  h«0.5ir.  In  the  figure,  l indicates  the  exit  length;  when 
i “0.01,  for  example,  it  indicates  that  the  leading  1%  of  the  foil 
has  exited  from  the  layer  of  water.  The  effect  of  the  thickness 
of  the  layer  on  the  pressure  distribution  during  exit  is  very  small. 
For  all  three  values  of  h considered  (h/tr*0.5,  1 and  2),  the  pres- 
sure distributions  are  very  similar;  when  X is  constant,  their 
differences  do  not  exceed  II.  Furthermore,  it  should  be  noted  that 
those  curves  shown  in  Fig.  10(c)  are  also  very  similar  if  both  the 
ordinate  and  abscissa  are  normalized  by  (l-£g).  In  this  case, 
only  three  curves,  one  for  each  value  of  X,  are  needed  to  describe 
the  pressure  characteristics  during  exit  of  the  circular-arc  foil. 
Such  a similarity,  however,  may  not  exist  for  other  foils  whose  gtt 
is  not  constant. 

4.  CONCLUDING  REMARKS 

Using  the  method  developed  in  [1],  the  problem  of  oblique 
entry  into  and  exit  from  a horizontal  layer  of  water  of  arbitrary 
thickness  by  a two-dimensional,  fully  ventilated  foil  has  been 
solved.  The  present  solution  contains  that  obtained  in  [1] . 
Therefore,  the  conclusions  derived  in  [1]  are  equally  valid  for  the 
present  case.  The  consideration  of  oblique  entry  makes  the  present 
solution  useful  in  studying  not  only  the  partially  submerged 
supercavitating  propeller  but  also  many  other  problems,  such 
as  the  high-speed  water  entry  and  exit  of  a wing  or  a control  sur- 
face, the  operation  of  a high-speed  hydrofoil  in  a high  sea  state. 
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Some  prominent  features  of  the  present  solution  are  summarized 
in  the  following. 

When  h is  infinite  and  the  flow  is  in  the  initial  entry  phase, 
analytical  solutions  in  series  forms  can  be  obtained  for  arbitrary 
foil  motions,  such  as  heaving  and  pitching.  In  other  situations, 
numerical  methods  are  more  efficient  to  effect  the  solutions. 

The  velocity  at  point  C (cf.  Fig.  1),  the  intersectional  point 
of  the  upper  free  surface  of  the  layer  and  the  cavity  wall  on  the 
suction  side  of  the  foil,  is  always  pointing  vertically  upward, 
irrespective  of  the  value  of  y so  long  as  0<y£*s.  The  magnitude 
of  this  velocity  is  equal  to  QQ(y)gt(0,0)cscYTr. 

When  gtt(s,t)  of  a foil  is  zero  during  the  exit  phase,  the 
present  linearized  theory  predicts  that  the  hydrodynamic  pressure 
acting  on  the  foil  is  zero  throughout  this  phase.  This  implies 
that  the  force  characteristics  of  an  exiting  foil  do  not  depend  on 
its  angle  of  attack  but  depend  on  its  camber  and  its  unsteady 
motions. 

From  the  numerical  results  obtained  in  this  work,  it  may  be  con- 
cluded that  the  effect  of  the  thickness  of  a layer  of  water  on  the 
force  characteristics  of  a foil  entering  the  layer  is  not  important 
during. both  the  initial  entry  phase  and  the  exit  phase.  However, 
the  thickness  of  the  layer  determines  the  magnitude  of  the  maximum 
loading  on  the  foil  as  it  goes  through  the  layer.  It  may  also  be 
concluded  that,  when  the  angle  of  attack  of  a foil  measured  from 
the  entering  path  remains  constant,  the  force  coefficient  -Cn  in- 
creases as  the  advance  ratio  X increases.  Nevertheless,  it  will 
never  exceed  the  force  coefficient  of  a corresponding  supercavita- 
ting  foil  at  zero  cavitation  number,  even  when  X ♦ » (that  is, 
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Figure  8(c).  Force  and  moment  coefficients  versus  irt/(hcscYi») 
for  a flat  plate  during  the  initial  entry  (0<t£l). 
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Distance  from  the  leading  edge  as  a fraction  of 

Figure  10(a).  Some  typical  pressure  distributions 
for  a flat  plate  foil. 
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Figure  10(b).  Some  typical  pressure  distributions 
for  a circular-arc  foil  with  6»a. 
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Distance  measured  from  the  leading  edge  along  the  chord 

10(c).  Pressure  distributions  on  a circular-arc  foil  (6=o) 
during  exit  from  a layer  of  water  of  h*0.57r. 


